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Abstract 


Errors  from  the  intended  performance  of  a  system, 
composed  of  specialized  co.  ponents,  could  be  due  to  any 
number  of  sub-system  fluctuations.  This  thesis  will 
examine  the  effect  of  these  resultant  errors  on  the 
mission  of  the  system  involved.  In  particular,  errors 
in  an  interceptor's  booster  performance  (seen  as  position 
and  velocity  errors  at  burnout)  will  be  propagated  along 
the  interceptor's  orbit  to  a  given  target  point  and  the 
resultant  error  at  the  intercept  point  will  be  found.  A 
range  of  intercept  orbits  will  be  examined  to  determine  if 
there  is  a  family  of  particular  intercept  orbits  that  give 
the  minimum  position  errors  at  the  intercept  point.  Con¬ 
clusions  will  be  made  as  to  what  typical  intercept  trajec¬ 
tory  should  be  flown  against  a  certain  type  of  target  so 
as  to  minimize  the  resultant  free  flight  error  propagation 
to  the  intercept  point. 
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I  Introduction 


Background 

In  the  problem  of  intercepting  any  orbital  target, 
there  are  always  inherent  abberations  in  the  interceptor 
system  components  which  cause  deviations  from  the  desired 
trajectory  during  the  course  of  the  intercept.  The  un¬ 
planned  errors  in  the  flight  path  could  be  due  to  any  number 
of  factors  which  could  effect  the  targeted  vehicle  in 
reaching  it's  programmed  orbit.  These  deviations  are 
caused  by  unavoidable  design  and  manufacturing  flaws  in 
the  components  that  compose  the  missle  booster  system,  such 
as  fluctuations  in  the  burn  of  the  rocket  motor  during 
orbital  insertion  or  small  errors  caused  by  a  slightly 
inaccurate  guidance  package.  All  of  these  factors  contri¬ 
bute  to  position  and  velocity  errors  which  occur  during  the 
powered  phase  of  the  mission,-^  -  up  to  booster  burnout. 
Hence,  the  errors  generated  during  the  powered  portion  of 
the  flight  will  be  the  ultimate  causes  of  errors  in  position 
and  velocity  of  the  interceptor  at  it's  planned  target 
point  termination. 
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Objectives 


This  study  will  be  concerned  with  investigating  how 
the  inherent  errors  of  a  missle  intercept  system  effect 
the  final  targeting  of  the  interceptor  involved.  Depending 
on  the  capability  of  the  booster  svstem,  each  interceptor 
has  the  versatility  to  reach  an  orbital  target  using  a 
variety  of  trajectories,  with  each  programmed  to  intercept 
the  body  at  a  different  point  in  it's  orbit.  Knowing  the 
accuracy  characteristics  of  the  booster  system  being  used 
allows  the  propagation  of  each  possible  intercept  trajectory, 
over  a  variety  of  times-of-f light ,  which  give  an  elliptical 
intercept  path.  Once  the  flight  path  is  traced  to  it's 
termination  at  the  target’s  orbit,  the  deviation  from  the 
proposed  intercept  point  can  be  found  by  comparing  the  ideal 
(no  error)  orbit  and  corresponding  intercept  point  to  the 
actual  orbit  found  through  propagation  of  the  error 
ellipsoid  from  booster  burnout  to  the  target’s  orbital 
path.  Thus,  given  an  error  covariance  matrix  modeling  the 
position  and  velocity  errors  at  burnout,  this  covariance 
matrix  can  be  propagated  by  means  of  a  state  transition 
matrix  to  the  predicted  intercept  point  to  give  a  covar¬ 
iance  at  intercept.  Then,  examination  of  the  eigenvalues 
and  eigenvectors  of  this  final  covariance  matrix  results 
in  the  generation  of  an  error  ellipsoid  for  the  interceptor 
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at  the  desired  target  point.  This  error  ellipsoid  can 
then  be  sliced  approoriatly  to  indicate  deviations  from 
the  intended  intercept  point  in  the  cross  track  directions 
of  the  target’s  orbit.  This  final  process  gives  actual 
miss  parameters  at  the  intercept  point  from  the  target  body 
in  question.  Repeating  this  process  over  the  range  of 
possible  intercept  trajectories  will  give  data  which  can 
be  examined  to  find  the  family  of  trajectories  for  which  the 
minimum  cross  track  intercept  errors  are  achieved. 

Research 

In  my  research  for  previous  work,  or  corresponding 
effort  on  my  topic,  I  found  no  data  or  documentation  on  the 
subject.  This  was  to  be  expected  due  to  the  uniqueness  of 
the  problem.  I  did,  however,  uncover  many  items  that  proved 
invaluable  in  the  development  of  this  project.  Several 
supportive  references  were  found  which  aided  in  the  extensive 
orbital  dynamics  modeling  which  was  done  in  the  program. 

The  typical  off-the-shelf  general  references  proved  invaluable 
in  helping  me  answer  the  "obvious"  questions  throughout  the 
project. 
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II  Computer  Simulation 


Assumptions 

A  computer  program  has  been  developed  to  simulate  the 
aspects  of  this  problem,  as  previously  outlined.  For  the 
purposes  of  this  study,  it  will  be  assumed  that  both  the 
orbital  target  and  the  interceptor  are  under  the  influence 
of  two-body  effects  only.  Since  the  typical  intercept  orbit 
would  probably  be  an  elliptical  one,  due  to  launch  and 
booster  constraints,  only  elliptical  intercept  orbits  will 
be  examined  to  further  simplify  the  process.  It  is  also 
assumed  that  the  satellite's  orbital  position  at  launch  is 
perfectly  known;  thus,  the  only  errors  present  at  the  inter¬ 
cept  point  will  be  due  to  position  and  velocity  deviations 
of  the  interceptor.  The  time  between  receiving  the  initial 
satellite  position  data  and  launch  is  assumed  to  be  zero, 
and  an  instantaneous  burn  time  of  the  interceptor’s  booster 
is  also  assumed.  In  order  to  simplify  the  implementation 
of  the  many  orbital  dynamics  formulas  needed,  a  universal 
variable  approach  will  be  used  throughout  the  project. 

Programming  Considerations 

Initially,  launch  site  data  will  be  inputed  to  find  the 
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position  of  the  launch  site  in  inertial  space.  Satellite 
data,  in  the  form  of  a  classical  element  set  for  the  target 
satellite  in  question,  will  also  be  inputed  to  find  the 
initial  position  and  velocity  of  the  satellite,  also  in 
inertial  space.  An  initial  time-of -flight  for  the  inter¬ 
ceptor  will  be  chosen  and  it's  corresponding  orbit  to  the 
intercept  point  (the  point  where  the  target  satellite  would 
be  after  that  time-of-flight  had  elapsed)  will  be  computed. 

A  test  will  be  made  to  see  if  this  computed  orbit  is 
elliptical.  If  this  orbit  is  not  elliptical,  then  the  time- 
of-flight  will  be  iterated,  and  the  initial  proceedure 
repeated  until  an  elliptical  intercept  orbit  is  found. 

Once  this  is  accomplished,  the  state  transition  matrix  for 
that  particular  elliptical  orbit  will  be  computed.  This 
state  transition  matrix  will  be  used  to  propagate  a  covar¬ 
iance  matrix,  representing  the  position  and  velocity  errors 
of  the  interceptor  at  the  booster's  burnout  point,  from  the 
burnout  point  to  the  given  intercept  point.  Thus,  a  covariance 
matrix  will  be  found  at  the  intercept  point,  representing  the 
propagated  position  and  velocity  errors  of  the  interceptor. 

By  examining  this  final  covariance  matrix,  conclusions  can  be 
made  as  to  the  effect  initial  booster  launch  errors  have  on 
the  final  miss  parameters,  associated  with  the  interception 
of  the  target  satellite. 


Computing  The  Launch  Site  Coordinates  In  Inertial  Space 


This  initial  phase  of  the  program  was  designed  to  allow 
the  earth  location  of  any  possible  site  to  be  entered,  as 
well  as  the  date  and  time  of  launch,  to  allow  flexibility 
in  the  final  simulation.  Inputed  data  consisted  of  the 
latitude  and  longitude  of  the  site,  it's  elevation  above 
sea  level,  and  the  year,  day,  and  time  of  launch.  A  sub¬ 
routine  was  developed  to  determine  the  local  sideral  time 
of  the  launch  site  given  the  site's  longitude  and  the  year, 
day,  and  time  of  launch  (Ref  ii 103-104).  Then,  knowing  the 
latitude,  altitude,  and  local  siderial  time  of  the  site,  a 
subroutine  was  developed  to  determine  the  position  of  the 
site  in  the  topocentric  reference  frame.  This  position 
vector  was  then  transformed  by  means  of  an  appropriate 
rotation  matrix  to  give  the  position  vector  of  the  launch 
site  in  the  inertial  frame  (see  Fig  I),  (Ref  It 98)* 

Computing  The  Satellite's  Initial  Position  In  Inertial  Space 

In  computing  the  satellite's  position  at  epoch,  or 
launch,  data  from  a  typical  element  set,  which  would  be 
received  at  the  launch  site  from  the  NORAD  space  track 
system,  was  used.  This  was  composed  of  the  length  of  the 
semi-major  axis  of  the  satellite's  orbit,  it's  orbital 
eccentricity,  it's  inclination,  the  longitude  of  the 
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Figure  1.  Relationship  Between  S-E-Z  and  I-J-K 


ascending  node,  the  argueraent  of  perapsis,  and  it's  true 
anomaly  ,  (see  Fig  2_).  Knowing  these  para¬ 

meters,  the  initial  position  and  velocity  of  the  satellite 
can  be  expressed  in  vector  form  in  the  perifocal  coordinate 
system  by  the  following  equations^/t^ 


where 


F  «  (r/  cos  V  P  +  fir*/  sin  V  <U 
p-  f^L  ~Sim  M  P  *  (g  +  C os  'U  )  ®.  J 

r>=  sCi-*.*-) 

in*  Vu  +  Cos  u) 


and  in  universal  variables 


»  X. 


Figure  2.  Orbital  Elements 


The  transformation  of  coordinates  between  the  perifocal 
system  (/’.O.V)  and  the  inertial  system  can  be 

accomplished  by  means  of  a  rotation  matix^(Ref  1«82). 
Thus,  if  are  the  components  of  a  vectorS.  in  the 

perifocal  frame,  then  the  coordinates  of  that  vector  in  the 
inertial  frame  are  given  byi 
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A  > 


•t 


where  > 


R,t  *  CoS  XL  c©s  CO  -  S' 0»o  jCL  9/a/  CJ  cos  X 

RiZ.  «  -  Co 3  JX  S/aj  Ci  -  S/aj  jO.  co 5  CJ  CoS  ^ 

/?  «3  ■>  5/*/  J/ai  xc 

tfz,  »  S/A/  jCZ.  Coy  6J  *  CoS  JX  J/At  6j  coj  X 

Rix  -  —  3//V  ^2L  57/v  CD  +•  Cos  JL  Cos  Co  CoS  X 

Re. S  *  ~  Co?  jCI  S/a/ 

R*t  -  _3//v  CD  S//e  X 

Rv*.9  CoS  (J  57*  ^ 

/?33  *  Coj  ^ 


Having  determined  the  elements  of  the  rotation  matrix,  it 
only  remains  to  find  f  and  P^in  terras  of  the  inertial  frame 

(y,^ 

inertial  frame  are  now  known  at  epoch  (Xa- launch  time). 


components  (x,3~,K).  Thus: 


r*  J  L/~*  1  L ^  J  Lv^J 


r 


<4N& 

Lr*  J  |>:  j 

Hence,  the  satellite’s  position  and  velocity  vectors  in  the 


Find  The  Satellite's  Position  At  The  Intercept  Point 

In  computing  the  orbit  needed  for  intercept,  the  time 
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of-flight  for  the  interceptor  was  chosen  as  a  variable 
which  would  give  a  resultant  elliptical  intercept  orbit. 

Due  to  the  fact  that  the  launch  time  was  chosen  to  coincide 
with  the  epoch  time  for  which  the  initial  position  of  the 
target  satellite  was  computed,  the  intercept  point  will  be 
the  point  on  the  satellite's  orbit  where  it  will  be  located 
after  traveling  the  chosen  time-of-flight  for  the  interceptor. 
Thus,  the  first  step  in  computing  the  intercept  orbit  is 
finding  the  intercept  point.  This  is  done  by  up-dating  the 
satellite's  position  from  it's  initial  location  at  epoch 
(■£=0)  to  it's  resultant  location  ,  marked  by  the  satellite's 
travel  during  the  time-of-flight  of  the  interceptor.  This 
prediction  problem  can  be  stated  as  follows* 


&  arc# : 

n 

0 

11 

£ 

Fi«o  : 

Jr 

,  V'  at  'rj/ye  ±r 

where  *£  is  the  time-of-flight  of  the  interceptor,  (see 
Pig  3). 

The  solution  consists  of  solving  the  Kepler  problem  for 
the  satellite  as  it  travels  from  K  to  P  in  the  given  time- 
of-flight.  The  initial  step  is  to  compute  the  universal 
variable  X  »  which  is  defined  by  the  relation  X-  Vr  . 
Kepler's  time-of-flight  equation  expressed  in  universal 
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Figure  3.  The  Kepler  Problem 

variables  is  soley  dependent  on  X  and  the  initial  position 
and  velocity  of  the  satellite,  as  seen  here* 

tqf~  +  (/-  £.)K*s*-  r.x  (7) 

where t  S  -  C  JT  -  JW  ST)/  IT* 

C  *  Cl-  Coj  ST)/  ? 

?  •  x  Vw 

An  intermediate  step  in  finding  the  radius  and  velocity 
vectors  at  a  future  time  is  to  find  \  when  time  is  known. 
Since  equation  (7)  is  trancendental  in  X  ,  a  trial  and 
error  solution  is  indicated.  Fortunately,  the  "6  vs 


il 


curve  is  well-behaved  and  a  Newton  iteration  technique  can 
be  used  to  solve  for  X  when  the  time-of-flight  is  given. 

If  a  trial  value  for  X  is  chosen,  call  it  X*»  then 

TOFn  *  (r\‘  +  o  -  $  +  r0XH  (8) 

where i  7b/*  is  the  time-of-flight  corresponding  to  the 
given  and  trial  value  X*. 

A  better  approximation  for  Xu  is  then  obtained  from  the 


Newton  iteration  algorithm* 


*--6, 


X«w  =  X*  IxsJCet 

.  (ft 

where*  t  is  the  given  time-of-flight  and  is  the 

slope  of  the  £  vs  X  curve  at  the  trial  point  Xn.  Now, 

knowing  X  we  wish  to  calculate  the  K  and  K  vectors  in 

terms  of  tZ }  VL  i  X  .  These  vectors  can  be  expressed  in 

terms  of  scalar  quantities  dependent  on  X  as  follows* 

r  *  -f  n  +  9  7.  (10> 

where*  _  .  ,  N  , 

fW-C^/rJc  £**-*’5 

■f xfes-/)/(ar)  J'J-x'c/r 

Thus,  the  satellite’s  position  and  velocity  up-dated  by  the 
chosen  time-of-flight  for  the  interceptor  is  now  known. 
Having  the  satellite’s  position  gives  us  the  data  for  the 
interceptor's  orbit*  the  initial  position  at  launch  and  the 
final  position  at  intercept. 


Figure  3«  The  Kepler  Problem 

variables  is  soley  dependent  on  X  and  the  initial  position 
and  velocity  of  the  satellite,  as  seen  here« 

TQF~  (/-  rmX  (7) 

where'  s-  (s*-,Mjr)/rr> 

C  =  (i-  TO")/  ? 

®  •  x  l/w 

An  intermediate  step  in  finding  the  radius  and  velocity 
vectors  at  a  future  time  is  to  find  X  when  time  is  known. 
Since  equation  (?)  is  trancendental  in  X  ,  a  trial  and 
error  solution  is  indicated.  Fortunately,  the  i  vs  )( 
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Compute  The  Orbit  Between  The  Launch  And  Intercept  Points 


At  this  point,  "both  the  launch  point  and  intercept 
point  of  the  interceptor  are  known.  It  is  now  necessary 
to  find  the  orbit  between  them.  In  other  words,  given  D  , 
(launch  site  radius  vector),  fit  (intercept  point  radius 
vector)  and  the  time-of-flight  from  Yt  to  fa,  ,  find  W, 
and  V*.  where  V*  is  the  velocity  vector  at  launch  needed  to 
intercept  at  Yl  and  Wi  is  the  interceptor's  velocity  at  fi  . 
Here,  the  time-of-flipiht  is  known  and  is  simply  the  chosen 
time-of-flight  for  the  interceptor.  Thus,  the  solution  of 
the  generalized  Gauss  problem  in  universal  variables  will 
give  the  needed  orbit.  A  subprogram  was  developed  to  solve 
the  Gauss  problem.  Initially,  from  Y.  and  Ifz.  ,  the  constant 
A  was  evaluated,  where s 

J >.  r-r  s/* 


A= 


(12) 


J  /-  CoS  A V 

Here,AV>  was  evaluated  by  finding  the  angle  between  and 
r*.  .  A  trial  value  for  2  was  then  chosen  and  the  functions 
S  and  C  were  evaluated  for  the  selected  trial  value  for 
2  where  the  series  representations  for  5  and  C  were  used 
to  eliminate  instability  problems  when  2.  is  near  zero. 


It 

1 

+!  6/ 

g3 

8  (  *  ... 

(13) 

J. 

?  4  2l 

* 

r/  7/  ” 

•  * 

"Z't  *  . .  . 

*  • 

(14) 

M 


.  1 


In  the  actual  program,  the  first  four  terms  in  these  series 
were  used. 


S  = 

c  = 


rr  -  j//v  jr"* 
/  —  CCS  J  ? 


Next,  the  auxiliary  variable  Y  was  found  from: 

V  y'-’i  4.  ^  0~*  ^ 

1  J“c 

and  ^  was  determined  from: 

Now,  the  trial  value  for  B  was  checked  by  computing  *£* 
from: 

■6.  =  x-5  S  *  A  5T~ 


(15) 

(16) 

(17) 

(18) 


(19) 


and  it  was  compared  to  the  desired  time-of-flight  which  is 
the  given  time-of-flight  chosen  for  the  interceptor.  If  it 
is  not  nearly  the  same,  the  trial  value  of  "fc.  is  obtained, 
within  some  convergence  criterion.  A  Newton  iteration 
scheme  for  the  adjustment  of  Z  was  used,  much  like  that  used 
//v  ~Hie  p*Bvicu%  section  f  and  was  implemented  as  follows: 


N+t 


where  "t  *  desired  time-of-flight 
is  given  in  equation  (19) 
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1 

1 


d*~ 

ZTk 


z:c  /4  8  C  c.  x  j 


S'-  it 

c'~  it  O-zwc) 


Here,  again,  instability  problems  could  occur  for  a  small 

value  of  Z  •  Thus,  series  expansions  for  5  and  C  were 

used  when  2  is  small  and  the  function  representations  for 

5*  and  Q.'  were  used  otherwise.  When  the  method  has  con- 

-  « 

verged  to  a  solution,  the  functions  ^  ,  and  3  were 


evaluated  from:  w 

f  =  /-  V. 

(23) 

ATT 

(24) 

Y/ri 

(25) 

Now,  since  fi-  •+  we  can  compute  V,  from: 

V>a  n--PH 

(26) 

and  similarly:  .  -*■  7* 

— »  S  •  • 

•  ^-5 — 

(2?) 

* 


Thus  now,  both  the  intial  position  (  )  and  velocity  (  V?  ), 

and  the  final  position  (  t\  )  and  velocity  (  Va.  )  of  the 
interceptor  are  known,  and  hence,  so  is  the  intercept  orbit, 
(see  Pig  4) 
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State  Transition  Matrix 


The  states  at  both  the  launch  point  and  the  intercept 

point  are  now  completely  defined  in  terms  of  inertial 

position  and  velocity  vectors.  The  transition  matrix 

between  these  two  states  can  now  be  found  by  taking  the 

partials  of  the  state  at  the  intercept  time  with  respect 

to  the  partials  of  the  state  at  the  intercept  time  with 

respect  to  the  partials  of  the  state  at  the  launch  time 

(Ref  3*129)«  n  ~  /.  \ 

c*  X  V,  c  )  i 

— ^  t  ^  <9  x  C"£  ‘  x. 


(z 


This  was  done  numerically  by  perturbing  each  state  element 
at  launch  by  a  small  amount  and  calculating  the  resulting 


new  orbit  over  the  same  time-of-flight .  The  difference 
between  the  final  perturbed  state  and  the  initially  known 
state  at  intercept,  divided  by  the  amount  of  the  initial 
perturbation,  gives  the  elements  of  the  state  transition 
matrix  between  the  two  states,  .  This  transition  matrix 
will  be  used  to  propagate  covariances  over  the  known  iner¬ 
tial  orbit. 

Covariance  Matrix  Propagation 

At  the  burnout  point  of  the  booster,  a  covariance 
matrix  was  developed  to  reflect  the  position  and  velocity 
errors  due  to  booster  performance.  These  error  parameters 
were  treated  as  input  to  the  program.  Thus,  the  position 
and  velocity  sigmas  were  inputed,  transformed  from  the  MKS 
to  the  universal  variable  system,  and  squared  to  give  the 
diagonal  elements  of  the  booster  covariance  matrix  at 
burnout,  which  is  itself  a  diagonal  6x6  matrix.  The 
upper  three  diagonal  elements  of  the  matrix  represented  the 
position  sigmas  while  the  lower  three  diagonal  elements 
represented  the  velocity  sigmas.  Then,  using  the  state 
transition  matrix  for  the  given  intercept  orbit,  this 
covariance  matrix  was  propagated  from  the  burnout  point 
to  the  intercept  point  along  tho  intercept  orbit.  This 
was  done  via  the  following  equation  (Ref  2)t 
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/ 


whcrs { 

P z  a  C oyiZjA/ucg-  AJsfTR/'K.  sir  ZLtoTern^epT 

P&  -  C  oi/s* A^A-rfH<  A--r  /3orno^j- 

Jji  -  Vc  /  R^tsS/  ^/G^  t7?A-rR/\ 

The  justification  for  the  use  of  the  above  transformation 


equation  is  seen  in  the  following  derivation  (Ref  3): 

0  f  .  _ _  ^  (pfFIN  lT'OH  \ 

o  X  '"^x)  ~  ft*,  ft  ;  \TkANnr  ion  M4trU) 


Px  =  £[3k(£x)SkG  Y]  ilo 


(TP/A/lT/oM  c>-f-  \ 

ce  /Wolfit  x  } 


0?) 


§o) 

6'> 


7~hus  o 

/x  -  2T JV  ft)  (Zftv0)  J*  ft)) 

"  ^  [*  (**>  *0  J *  ft)  i*  ftrx  ft  ft)TJ 
-  Ufa,**)  E [jvfrx)  J <  (^)T]S(%-i.)T 
Pr  ~  JT  ftft)  P3  jITftft )*7'  (32) 


where «  /T£  3  indicates  the  expected  value  of  the 

braclcetted  quantity.  Since  this  thesis  is  interested  in 
only  the  resultant  position  errors  at  intercept,  the  pro¬ 
pagated  covariance  matrix  was  examined  and  the  upper  3x3 
sub-matrix  was  extracted.  This  sub-matrix  represents  the 
propagated  position  errors  at  intercept  since  the  initial 
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4 


/ 


covariance  matrix  was  developed  by  placing  the  position 
error  elements  in  the  upper  diagonal  locations  of  that 
matrix.  This  position  sub-matrix  represents  the  position 
errors  at  intercept  with  respect  to  the  interceptor's 
orbit?  «  the  deviations  from  the  intended  intercept 
point  referenced  from  the  track  of  the  interceptor's 
orbital  path.  However,  the  items  of  interest  are  the 
deviations  from  the  location  of  the  intended  target 
relative  to  that  target.  Thus,  the  position  sub-matrix 
requires  a  transformation  so  as  to  reflect  the  actual 
position  errors  referenced  to  the  relative  velocity  vector. 
This  was  accomplished  by  a  rotation  of  the  position  covar¬ 
iance  matrix,  rotating  about  the  relative  velocity  vector 
between  the  satellite  and  the  interceptor.  This  was  done 
by  realigning  the  basis  vectors  of  the  error  ellipsoid 
along  the  relative  velocity  vector  between  the  inter¬ 
ceptor  and  target  satellite  (see  Fig  5).  This  rotation 
matrix  was  developed  as  follows  (Ref  li6l)« 


3  ® 
StN 


C05  ^ 

4 


Situ  &  Cos  4 


Cos  Q  $IN 

Cos  4 

% 

S/tt  Q 


-  S/n  e 
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Figure  5*  Slice  Ellipsoid  To  Obtain  Relative  Error  Axes 


Using  this  rotation  matrix,  the  position  covariance  matrix 
was  rotated  into  the  relative  velocity  vector  as  follows 
(Ref  2 ) «  Given  initial  covariance  elements  and  desired 
rotated  covariance  elements  Jx'  i 

p'=  p*  e[s*  jvj 

p‘=  elR's <  CK'Ji^T] 

=  e  [r'jz 

=  P'e[j*£-K  TJ«'T 

p'  =  P  R,r  Cr?) 


Now  the  position  error  ellipsoid  at  the  intercept  point 
relative  to  the  satellite  is  known  (see  Fig  6).  It  now 
remains  to  find  the  actual  miss  distances  from  this  covar¬ 
iance  matrix.  Due  to  the  intercept  problem  involved,  the 
in-track  errors,  those  errors  that  indicate  if  the  inter¬ 
ceptor  is  either  in  front  of,  or  behind  the  target  at 
intercept,  are  of  no  real  interest  to  this  study,  for 
they  simply  dictate  if  the  intercept  should  occur  a  little 
sooner,  or  later,  than  expected.  Thus,  the  row  and  column 
of  the  position  covariance  matrix  that  reflect  the  in-track 
errors  were  eliminated,  leaving  a  2  x  2  matrix  representing 
the  cross-track  errors  at  intercept.  The  eigenvalues  of 
this  matrix  were  then  found  and  their  square  roots  were 
taken  to  determine  the  axis  lengths  of  the  cross-track 
error  ellipsoid  at  the  intercept  point.  Once  these 
principal  axis  lengths  were  converted  back  to  the  MKS 
system,  these  lengths  represented  the  propagated  miss 
distances  from  the  target  point  (in  meters).  These  error 
ellipsoid  axis  lengths,  representing  the  error  from  the 
center  of  the  ellipse  to  a  point  on  it's  edge,  are  one 
sigma  errors.  They  represent  a  68%  confidence  level. 

Dufe  to  the  linearity  of  the  error  propagation  process, 
the  three  sigma  or  99 %  confidence  level  can  be  found  by 
simply  multiplying  the  errors  by  a  factor  of  three. 


a 


Thus,  the  error  ellipsoid  footprint  cross-track  to  the 
relative  velocity  vector  is  now  known,  along  with  the 
magnitudes  of  each  perpendicular  arm  of  the  axis  which 
make  up  that  cross-track  ellipsoid.  Throughout  the  deve¬ 
lopment  of  the  program,  the  reliability  of  each  major  sub¬ 
program  or  element  of  the  main  program  was  examined  in  a 
variety  of  ways.  All  of  the  program  elements  concerned 
with  orbit  determination  were  verified  by  using  data  from 
examples  previously  done  in  the  references  from  which  that 
element  was  developed.  Only  high  correlation  between  the 
output  of  that  sub-program  and  the  “book  answer"  insured  the 
proper  functioning  of  that  program  element.  The  matrix 
manipulations  in  the  latter  part  of  the  main  program 
were  verified  by  examining  output  matricies  during  inter¬ 
mediate  steps  and  calculating,  by  hand,  certain  elements 
of  the  matrix  found  in  the  subsequent  step  and  comparing  the 
corresponding  numbers.  This  technique  assured  the  reliability 
of  the  intricate  matrix  manipulations  used  in  the  propagation 
of  the  initial  covariance  matrix  to  the  intercept  point. 


Ill  Program  Execution 


With  the  theory  now  modeled,  the  program  is  now  ready 
to  produce  the  data  required  from  which  appropriate  con¬ 
clusions  can  be  made.  The  launch  site  was  held  constant, 
and  a  variety  of  target  orbits  were  examined  to  make  up 
one  senario.  For  each  orbit  examined,  varrying  times-of- 
flight  for  the  interceptor  were  used  until  an  elliptical 
family  of  intercept  trajectories  were  found.  Then  the 
time-of-flight  was  slowly  iterated  over  this  elliptical 
range  to  produce  a  quantity  of  intercept  points  and,  hence, 
two-dimensional  position  error  data  at  each  point.  This 
array  of  data  was  examined  and  the  minimum  horizontal  and 
vertical  position  errors  were  extracted  along  with  their 
corresponding  time-of-flight,  which  indicates  on  what  type 
of  intercept  trajectory  they  occured.  Also,  plots  of  the 
cross-track  error  ellipsoid  foot  prints  were  made  for  the 
first  and  last  elliptical  orbits  in  the  family  and  also  for 
the  cases  of  the  minimum  horizontal  and  vertical  miss 


distances. 


Results 


In  implementing  the  program  to  generate  data,  many 
target  orbits  were  examined.  For  all  the  runs  used,  the 
launch  site  was  fixed  along  the  equater  with  zero  longitude. 
This  located  the  site  along  the  first  point  of  Aires,  T  . 

long  the  X  unit  vector  in  inertial  space.  Due  to  the 
restriction  on  the  geometry  of  the  intercept  orbit  (it  must 
be  elliptical),  many  of  the  target  orbits  considered 
yielded  too  little  information  to  be  examined  for  definite 
conclusions.  This  was  done  to  the  fact  that  there  were  very 
few  elliptical  intercept  paths  that  could  be  flown  against 
them.  Thus,  in  this  section  only  two  types  of  target  orbits 
were  considered.  A  direct  polar  target  orbit  which  passes 
perpendicular  to  the  site  was  used  along  with  a  moderatly 
inclined  orbit  which  also  passes  perpendicular  to  the  site 
(see  Fig  7),  (see  Table  I).  Initially,  cross-track  error 
data  was  obtained  by  varrying  time-of-flight  for  the  inter¬ 
ceptor  and  up-dating  the  satellite's  position  from  epoch  by 
that  time-of-flight.  This  technique  gave  the  range  of 
elliptical  intercept  orbits  over  a  segment  of  the  satellite's 
orbital  path,  the  eccentricity  of  each  target  orbit  was 
set  at  zero  (circular  orbit)  and  .2.  (moderatly  eccentric) 
to  obtain  data  sets  for  comparison.  Initial  covariance 
values  at  burnout  were  found  by  inputting  sigmas  of  10  meter 
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Table  I>  Target  Orbits  Examined. 


Target  Orbit  #  1  ••  Target  Orbit  #  2: 


» /.  out,  IS 

1.07 /t  is 

.  /  <  /0“fi 

<D.Z 

i  ?o° 

A.-  ?0° 

=.  270° 

JTL  «  270* 

-  £  “ 

CO  *  ©* 

*  O  ° 

V  *  0° 

:et  Orbit  #  3« 

Target  Orbit  #  4) 

«  l onus' 

r  /  x  /ST* 

C  65-° 

* 

*  o  * 

-  Z7oQ 


*  /  07 /6/5~ 
<9.Z 
Z*  65"" 

.a  «-  9a* 

63  -  fi>* 

M  •-  Z7©* 
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in  position  error  and  . 030<f3  meters  per  second  (  ./  foot 
per  second)  in  velocity  error.  Then  two  other  data  sets 
were  developed  by  first  zeroing  out  the  position  error  at 
burnout  and  then  the  velocity  error.  This  was  done  to 
examine  the  sensitivity  of  the  resultant  intercept  error  to 
burnout  position  and  velocity  errors.  This  cross-track  error 
data  for  each  target  orbit  is  summarized  in  table  II. 

Plots  of  this  data,  cross-track  error  vs.  tine-of-flight , 
can  be  seen  in  Figures  (8  to  24).  The  parameter  "circular¬ 
ized  cross-track  error"  was  developed  to  combine  the  two 
axis  errors,  as  the  other  plots  show  both  error  ellipsoid 
axes  plotted  as  they  propagated  vs.  time.  It  was  found  by 
determining  the  radius  of  a  circle  with  the  same  area  as 
the  ellipse  represented  by  the  two  axis  lengths  of  the 
cross-track  slice  of  the  error  ellipsoid  at  the  intercept 
point.  Next,  identical  data  computation  scheme  was  used 
to  examine  error  propagation  when  just  one  point  on  the 
target  orbit  was  used  as  the  intercept  point.  Here,  the 
satellite's  position  at  epoch  was  determined  and  held  con¬ 
stant  as  a  family  of  elliptical  intercept  orbits  were  com¬ 
puted  for  various  times-of-f light .  These  intercept  orbits 
were  used  in  the  error  determination  to  examine  the  effect 
of  orbit  geometry  on  the  intercept  error  growth.  The 
corresponding  data  range  is  shown  in  table  III. 
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Table  II i  Intercept  Errors  Over  Range  of  Target  Positions. 


Target  #  1. 

Booster  Burnout  Errors i  Position  /O  M 

Velocity  (P.C>30$8 

Variance  in  the  Principal  Axes  Lengths  (M)*  2.1- "*3  4$, 5S- 

4S.OI  -  /2o.  S9 

Variance  in  the  Circular  Error  Parameter  (M)«  3/. 27 -74.  92 

Booster  Burnout  Errors*  Position  /x/0“*  M 

Velocity  0  0  ^/scc. 

Variance  in  the  Principal  Axes  Lengths  (M)s  ^ 

39.  //  -  //Z.84 

Variance  in  the  Circular  Error  Parameter  (M):  23.  7/.  89 


Booster  Burnout  Errors*  Position  / O  *1 

Velocity  /*  /0~*  “fee 

«4. Of"  3.88 

Variance  in  the  Principal  Axes  Lengths  (M)* 

2.2.39  -  43. 

Variance  in  the  Circular  Error  Parameter  (M)*  9 .37-/1.44 


Target  #  2. 

Booster  Burnout  Errors* 


Position  /O  /» 
Velocity  0-0  3048 


Variance  in  the  Principal  Axes  Lengths  (M)» 


20.35*  -  48.4fc 
39.  /7  -  /0  7.  9Z. 


Variance  in  the  Circular  Error  Parameter  (M)»  2,3.23  -  67.65' 
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Booster  Burnout  Errors: 


Position  / * /O  ^  A7 
Velocity  0>0  3098  '’Vjfcc. 
Variance  in  the  Principal  Axes  Lengths  (M): 

za.23-  47.65- 

Variance  in  the  Circular  Error  Parameter  (M):  2.if.  79  -  6?.  7/ 

Booster  Burnout  Errors:  Position  /O  Af 

Velocity  /*  /o“8  ^/sec 

.  .  4.76-  3.09 

Variance  in  the  Principal  Axes  Lengths  (M):  __  „ 

/9.77  ~  3  /.  O/ 

Variance  in  the  Circular  Error  Parameter  (M):  9.7^-  9.74 


Target  #  3. 

Booster  Burnout  Errors : 


/<?  A* 

0.0  309  8  /*^C 

2.0.72  -  90.96 
2.6.43  -  /52. 89 

Variance  in  the  Circular  Error  Parameter  (M):  23.4/ —78. 63 


Position 
Velocity 

Variance  in  the  Principal  Axes  Lengths  (M): 


Booster  Burnout  Errors! 


Position 

Velocity 


!*/ o~s  r* 


o.oso4g  'H/sexi 
,  s  /9.9V-  3?.  76 

Variance  in  the  Principal  Axes  Lengths  (M):  ^  c/  _  /40  SB 

Varinace  in  the  Circular  Error  Parameter  (M)«  z/.  /&  -  79.  8f 


Table  II  (cont.) 


Booster  Burnout  Errors:  Position  /0 

Velocity  /*/0‘*  y*1/s'<fc 

6.44  -  4.4/ 

Variance  in  the  Principal  Axes  Lengths  (M)i  ^  ZQ  -  f?  S7 
Variance  in  the  Circular  Error  Parameter  (M):  7.2  4-/6.22. 

Target  #  4. 

Booster  Burnout  Errors:  Position  /Om 

Velocity  0.030*8  '"/sec. 

Variance  in  the  Principal  Axes  Lengths  (M):  ^'4 2-  ‘/SQ  08 

Variance  in  the  Circular  Error  Parameter  (M):  *74:,  <£9 

Booster  Burnout  Errors:  Position  /  «  /*7 

Velocity  0-030*8 

/7.  Zt  -  34.  yz 

Varinace  in  the  Principal  Axes  Lengths  (M):  7^-,  /Z7.  o*? 

Variance  in  the  Circular  Error  Parameter  (M):  /8.84-  70.87 

Booster  Burnout  Errors:  Position  /0-^f 

Velocity  /<< /O  8 

6.7s--  f.zz 

Variance  in  the  Principal  Axes  Lengths  (M):  ^ 

Variance  in  the  Circular  Error  Parameter  (M):  9.43-  IS  IS 
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Table  #  III:  Intercept  Errors  for  a  Single 
Family  of  Interc ept  Orbits. 


Booster  Burnout 

Variance  in  the 
Variance  in  the 

Booster  Burnout 

Variance  in  the 
Variance  in  the 

Booster  Burnout 

Variance  in  the 
Variance  in  the 


Booster  Burnout 

Variance  in  the 
Variance  in  the 


Target  #  1. 

Errors:  Position  /Q 

Velocity  <D'0"$o4 & 
Principal  Axes  Lengths  (M): 


2.2.03  -3 n.  ys- 
**4.  *7  -  /FS-.17 
Circular  Error  Parameter  (M):  -  7g 


Error: 


Position 

Velocity 


/  *  /o'8 

o.eio+% 


Principal  Axes  Lengths  (M):  38.86 

3fc.46~  /44.07 

Circular  Error  Parameter  (M):  Z^.7z  - 


Errors:  Position  lO 

Velocity  /*  /o"8 

^•27  "^.6/ 

Principal  Axes  Lengths  (M):  2/  _  yy 

Circular  Error  Parameter  (M):  9.64 -/6.3V 


Target  #  2. 

Errors:  Position  / ©  -'+» 

Velocity  0  ^/sco. 

„  .  .  ,  .  T  ...  /e.4-2  -  35*.// 

Principal  Axes  Lengths  CM) *  _ ^ 

35".  "74  -  /JY.  *  8 

Circular  Error  Parameter  (M):  £5-  ^  <^g 
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Table  #  IIIi  (cont.) 


Booster  Burnout 

Variance  in  the 
Variance  in  the 

Booster  Burnout 

Variance  in  the 
Variance  in  the 


Booster  Burnout 

Variance  in  the 
Variance  in  the 

Booster  Burnout 

Variance  in  the 
Variance  in  the 


Errors:  Position  /  k  /O  sh 

Velocity  0.©3o<*B  '"/sec. 

Principal  Axes  Lengths  (M)»  *7.y4-  34.5'f 

3 /.  f9  -  n *.  *7 

Circular  Error  Parameter  (M):  23.07-  69.38 


Errors:  Position  /€)  svi 


Velocity  /*  /0~8  '"/sec 


Principal  Axes  Lengths  (M): 


f.ii  -  3.98 
/7.00  -  5^.83 


Circular  Error  Parameter  (M):  /j-.fZ 


Target  #_3. 


Errors:  Position 


IO 


Velocity  0-0 
Principal  Axes  Lengths  (M) 


304-8  ^/sec. 

/9..T9-  37.93 
8  23 .*8-/4o.*Z 


Circular  Error  Parameter  (M):  2/.<ty  -  72.98 


Errors:  Position  /x  /o-^* 

Velocity  ©.o3©*8  "**/ Sec 

Principal  Axes  Lengths  (M):  /S-/3  -  37.62. 

70  /3  -  '29.  62. 

Circular  Error  Parameter  (M):  /7./0- 69.^7 
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Table  #  III»  (cont. ) 


Booster  Burnout 

Variance  in  the 
Variance  in  the 


Booster  Burnout 

Variance  in  the 
Variance  in  the 

Booster  Burnout 

Variance  in  the 
Variance  in  the 

Booster  Burnout 

Variance  in  the 
Variance  in  the 


Errors  i 


tC  * 

/  *  /o‘*  '"/sec. 

*7.25*-  2 .8* 

/Z.Zc-  sr.ss 

Circular  Error  Parameter  (M):  ^.«?©  - 


Position 
Velocity 
Principal  Axes  Lengths  (M)» 


Target  #  4. 

Errors:  Position 

Velocity 
Principal  Axes  Lengths  (M): 


/O  St* 

o. ©3048  '"/sec. 

/fc.57-  36.23 
2*. 42-  /37 .6/ 


Circular  Error  Parameter  (M):  £2.£>a-  73. Z 9 


Errorsi  Position  /  *  /0“% 

Velocity  o 

.  .  -  36. oy 

Principal  Axes  Lengths  (M)«  „  ,  „ 

2.®. 72-  /Zt.?7 

Circular  Error  Parameter  (M):  /<?,  3r- ^2. 37- 


Errors:  Position  IO 
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FIGURES  12-15. 


Data  for  the  Range  of  Target  Positions 


TARGET  #  2. 
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FIGURES  16  -  XI- 
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CROSS  TR^CK  ERROR  (M) 
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Figure  #  19 
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CROSS  TRRC  <  ERROR  (M) 
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FIGURES  24-2 
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CROSS  TRRCK  ERROR  ( M ) 
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CIRCULARIZED  CROSS-TRACK  ERROR  CM3 
J?4.00  92.00  40.00  48.00  56.00  64.00  72.00  80.00 
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Figue  #  31 
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FIGURES  32  -  35. 


Single  Family  of  Intercepts 


TARGET  #  3. 
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CROSS  TRACK  ERROR  ( H ) 
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FIGURES  36  -  39. 


Single  Family  of  Intercepts 


TARGET  #  4. 
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TROSS  TRAC  <  ERROR  ■  M ) 

_f\Or  *O.OC  4CK  or  CO. 00  80.00  100.00  120.00-  140.00 

■  *  i  i  •  i  i  i 


Figure  #  j6 
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Conclusions 


In  analyzing  the  results,  the  first  obvious  con¬ 
clusion  is  the  fact  that  shorter  times-of-flight  give  the 
smallest  miss  distances  at  the  intercept  point.  This 
observation  is  consistant  throughout  the  data  analysis 
and  does  not  depend  on  orbit  geometry.  But  this  fact 
was  an  expected  result  in  that,  due  to  the  linearity  of 
the  propagation  routine  and  a  bit  of  common  sense,  the 
errors  should  increase  as  the  interceptor  tracks  further 
down  it's  orbital  path,  as  the  times-of-flight  become 
larger.  In  comparing  like  target  orbits  which  have  diff¬ 
erent  eccentricities  it  was  seen  that  the  eccentric  orbit 
yielded  smaller  intercept  errors  for  the  same  time-of- flight. 
The  length  discrepencies  appeared  in  both  principal  axes. 

This  fact  indicates  that  the  resultant  intercept  error  is 
dependent  on  the  eccentricity  of  the  target  satellite's 
orbit.  Comparing  data  for  identical  times-of-flight 
between  the  two  different  geometry  orbits  considered, 
it  was  evident  that  the  inclined  target  orbit  gave  the 
smaller  intercept  errors.  This  suggests  that  the  inclina¬ 
tion  of  the  target  satellite  has  a  definite  effect  on  the 
final  target  miss  distance.  But  a  much  more  extensive 
examination  along  this  line  should  be  attempted  before  a 
definitive  statement  concerning  this  effect  can  be  made. 
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The  most  generalized  observation  that  can  be  made  in 
examining  the  results  concerns  the  magnitudes  of  the 
computed  errors  themselves.  The  fact  that  these  one 
sigma  errors,  represented  by  the  circularized  error  para¬ 
meter,  typically  range  between  2-/3  and  7<?:9meters  indicates 
that  the  final  free  flight  intercept  position  is  known 
better  than  the  position  of  the  target  itself.  For  the 
purposes  of  this  study,  the  target's  position  at  intercept 
was  assumed  to  be  perfectly  known  when  in  reality  it's 
three  sigma  location  is  only  known  to  the  order  of  about 
500  meters.  Now  comparing  the  same  target  orbit  data 
while  initial  position  and  velocity  errors  were  nulled 
shows  that  the  final  intercept  error  is  mainly  dependent 
on  initial  booster  velocity  errors  and  is  not  very  sensitive 
to  corresponding  position  errors  at  the  burnout  point. 
Finally,  in  examining  the  plot  of  the  circularized  error 
parameter  vs.  time-of-f light  it  is  noticed  that,  for  all  the 
cases  considered,  it  is  a  straight  line  increasing  over  time. 
This  result  was  also  expected,  as  stated  previously,  due  to 
the  linearity  of  the  error  propagation  over  time  and  this 
fact  lends  considerable  credibility  to  the  computed  data. 


Recommendations 


Throughout  the  development  of  this  thesis,  and 
especially  in  the  final  stages,  I  became  more  and  more 
aware  of  the  many  things  that  could,  and  should  be  done, 
for  which  I  just  would  have  time.  The  next  step  in  the 
utilization  of  the  program  would  be  to  alter  the  logic  so 
that  the  covariance  propagation  is  done  from  an  actual 
burnout  point  to  the  intercept  point.  Because  of  the 
initial  assumption  of  an  instantaneous  burn  of  the  booster, 
the  program  propagated  the  covariance  from  the  launch  site 
to  the  intercept  point.  An  appropriate  burnout  altitude 
could  be  inputted  as  a  variable  and  it's  location  on  the 
path  of  the  intercept  orbit  could  be  found.  It  then  remains 
to  use  this  location  to  begin  the  covariance  matrix  pro¬ 
pagation  to  give  more  realistic  results.  Another  interesting 
area  which  could  be  examined  is  the  development  of  a  contour 
map  of  time  of  launch  for  the  interceptor  vs.  cross-track 
error  for  various  t imes-of-f light .  Considering  a  certain 
launch  time,  corresponding  launch  and  intercept  points 
could  be  found  and  a  family  of  intercept  orbits  computed 
between  the  points  along  with  the  associated  miss  distance. 
Then  the  launch  time  could  be  increased  and  the  process 
repeated  until  a  wide  range  of  launch  times  have  been 
considered.  Appropriate  conclusions  could  then  be  made  as 


to  the  effect  of  the  time-of-  launch  on  the  ultimate  inter¬ 
ceptor's  error  at  the  target  point.  Also  of  interest  might 
be  the  examination  of  the  sensitivity  of  the  interceptor's 
error  to  intercepts  on  different  passes  of  the  target  rel¬ 
ative  to  the  launch  point.  Here,  hopefully  conclusions  could 
be  made  as  to  which  revolution  of  a  target  should  be  used 
to  fly  the  most  accurate  intercept  trajectory  against. 

In  the  implementation  of  this  program,  I  see  that  I've 
just  scratched  the  surface  by  finding  the  cross-track  errors 
at  the  intercept  point  and  doing  the  rather  general  orbit 
geometry  comparison.  This  groundwork  has  been  laid  in  hopes 
that  follow-on  research  would  be  attempted  to  further  broad¬ 
en  the  knowledge  on  this  subject. 
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APPENDIX  A 

COMPUTER  PROGRAM 


79 


o  o  h-  no 


5|«0  GRAM  SAT  (IN|OliT  =  /«0,0UT:sjr,T  APES = OUTPUT ) 

DIMENSION  ‘'(3,6)  ,e(E,fc)  ,C(5,i>  ,0(3,1)  ,SS.'3,1)  ,WK(60>  , CHECK  (0,6) 
DIMENSION  ®HI (6,6)  ,RR(3,3) , ROT  (3,3) , COV (6,6) ,C0VR(2,2) , DD(3,3) 
HT PENSION  Er ( 3 , 3 )  , (I (2, 2), 4(2) 

01  PENSION  VAT  Pf>  ,ZA(400>  ,ART  (810)  ,TT<4  90)  ,  OEP (4 0 0) 

COMMON  THETA 

COMMON/ LSI  IM-/LPNGI,DAV,Tr  “IE,  THETA  3D 

COMMON/ sitee/lati, HI 

CCMMON/TIKE/T 
COHPL  EX  1 J  , H 

REAL  LATT,lONGI,LAT,LONG 
*lAL  NJ,MtGl,MAG2,Ll,L2 
INTEGER  CP 


0  *  4  444*444  *4*141  4**X*-***X*4*4*,*»*4  4*»  **,***4*1 
04*4*4  4  **»,-***«  ****.*•» 

C*v44  ERROR  PROPAGATION  FDR  GROUND -TO- SATELLITE  ELLIPTICAL  *• 

c.,*44  INTERCEPT  QROITF  JO  EXAMINE  O.E.p.  GROWtH  OF  BOOSTER  COVARIANCE  ** 

0**4  4.  4  *4-  4X4**X  *«  4*4**. >***4**4*4*  ,**4***»  4,4** 
^4*444+*4X»*».**»**4****»******  +  **44»***44***4«**J 


‘***  ******** 


PRINT  98 
PRINT  38  - 

8b  P  ORMAT ( IX , 1 23 ( IP* ) ) 


11 


KLM=0 
1  JI  =  G 

IF (III. EO. 3) GO  TO  222 
KLM=KLM+1 


111=111+1 
C0  =  0 


LK=CI 

PRINT*,"  ” 

PRINT*,  ORBIT  NUHOER  ",III,"  (MMM )##** 

PRINT* , "  ” 

4*4*4*. »**»»*,x,44***4********4,*,*4**4,  4*4**  *4* 


*  *♦*  4***4 


INPUT  OATt 


0*41  »»»»»*»(.>  »»******»»< 

REAO*,T,IF, JF 
READ*  ,RP,  \l° 

READ*  ,L  ATT,  L3  NO’  ,HI,  TIME,  3W,  THETAS  0 


n..*.****-****#*.*!*  •***»*•* 


4»4f.»4.*,.  **.»*-  l*4,i 


PRINT*,  "COVARIANCE  VALUE1?..  .POSITION  £PROR=" , RP , " < 1 >" 

print*,"  velocity  erpor=",vp," (m/seo 

PRINT** "  •• 

RP=(RP/6  r  81^!?)  **2 
VP=(VP/79f:?.366?8)*'‘2 
00  15  1*1,3 
IE-  C0V(I.I)=P3 

no  16  J  =  4,5 
lfc  COV(J.J) =v° 

no  9  1=1,  A  00 
YA  (I)=0. 

1  A(I)=0. 

ARY(I)=0.  , 

9  TT ( I) =3  . 

****»**1>1*44»44>*4>4  4  .  .  4*4  A**.*******.**.**** 

C  COMPUTE  LOCAL  SJDERIAL  THE  OF  LAUNCH  SITE 

014  1*1  *444*44441  4*444>444>444444*4f4*44  4>Mt*>  4*44444144414444 


M  1  4*4 


4444l>4«4*> *»' 


‘4441444 


CALL  LSTIM" 


‘•444444' 


4  4  4  444  ‘  ■ 


T=  .1 
THETA=0 . 

04*l«,4>*X4444tl*44r.«x*44*.44444  44444444  4  44  4  444*441 

C  COMPUTE  INERTIAL  COORDINATES  OF  _  A  JNCH  SITE 

0**l4lt*4*<44*4|  4»4»4.r*44.4*»*4^^.'*f4*»*l*44»*****4**  4444 

CALL  SITE(X ,Y , 7  ,  VX  ,  V Y  ,  VT) 

0*4*4**444444441  144 4*4 4«4< 444**4 44  *  4 44 4  4  4  »  * 44» 444  I 

C  COMPUTE  INITIAL  SATELLITE  POSITION 

Q*  •  *♦♦•».**♦*  *14  *4  14*  *4*  4*1  ♦♦♦**  •**  4441  *  *«  *  4  44  *♦  »  4  ***♦*»♦*  *1 

CALL  ORDIMXI,  Y  T  ,  ZI ,  V  Y I ,  Vf  I  ,VZI,RO,  VO,  A) 

u=a 

j  jj=i 


1M4411M4**'  *' 


>41*44  441  lt>  ***>41 


.  4  4444444*444 


o  o 


UPDATE  SArELI_rTE  POSITION  10  INTERCEPT  c0INT 

£*».».  »  .  ♦•»*  4  »  ♦  *  »  V  f  Mf 

D0T=XI*V*  I*YI*  VVI  +  7J*  v/7I 
7 1  CALL  FINny<XX,nOTfR3tVOtIO>*> 

IFCUD.  ED. 50)  DO  TO  19 

CALL  LOCArE(XX,PO,VO,F,G,S*C:,Z?, *> 

Rl  =  CfcXI  *G'  \/  XI 
R2  =  FJ,YI*0*  VYT 
R3  =  F+  71  +Gr  VTT 
0c=xx*(7t.5_1> ) /(R*?o) 

'JG  =  1 .  -XX**  2 ,‘CC/R 

V<-  =  0F*XI*PG,'\/XI 

VT  =  0F*YI+PS*vfYI 

VC=DF'7I+D3 -77J 

R=S0RT<K1» *  2. +  F  2**2. *P3**2.> 

V=S0RT( V3*  ‘  2*  VG-’  *  2+V4**2> 

-r  ■»  4*4  .4*4****<  *  *  *  X  *  4  4»4*****MHU»»H  *4  *  *  *  '■  44  4  4 

C  COMPUTE  OEDIT  PfOUIRED  FOR  =>RDRE  INTERCEPT 

MAG1=S0RT ( X*  *  2 • *Y**2. ♦?**?.) 

M  A  G2  =  SDRT  ( R 1*  *  2  .  «-R2*  ♦  2  .  *R3*  *  2  . ) 

Xl=X/MAGl 
Y1=Y/YAG1 
Zi=7/MAG1 
X  2=R1 /MAG2 
Y2=R2/MAG? 

72  =  R3/-1  AG2 

n0P=Xl»X?4Vi* Y2+Zl*72 
NU=ACOS (CD3) 

CALL  GAUSS(M'J»MAGl*MAG2»r»GfDG»L^) 

IF  (LM, £0.1) GO  T0  77 
IF  (LM.EO.f ) GO  TO  77 
'Jx  1=  (pl“Fr  X )  'G 
7I2=(R2-F* Y) 75 
VI  3=  (R3-c‘  7  )  7 G 


7T1  =  (o;;»ri-x>  fr, 

VT2=(0G*R?-Y) /G 
\/7  3=  <  DG*  ^3-^) /G 

vi=sopt  (vnM’ .  +vi2**2.  +'/r?*  *2 .) 

77S=SQR7(  7T1**  2+7T2»*  2  +  VT  3 *  *2) 

7VI=SQRT  (7Il*‘2.  +  7I2**2.«-7I3+*2.,> 

0071=X* 7I1+Y* 71 2+Z*Vl3 
CC=7I+*2.-l .  /HAr-i 
E1=CC~X-DDT1*VI1 

E2=CC*f -DOT 1*712  - 

E  3=CC’ 7-0DT1*  V 1 3  f 

EIJ=SQ*n  (El**?.  +E2**2. +E3**2.) 

71=7  .  1 

4M  =  Af’S(l./(VI**2-2./MA6in 

C'**1'*'**-'  ****44**«  **♦♦♦<  fc**.*.*.*****.***.*.***-*  4»44»»4  44‘*'»44U| 

C  OPANICH--ELLI;5rirAL  OR  HY=E?3DLIC  INI  ERCEDT  OkOIT 

0  •»  4  4-*  +  M*«4**t  4  »  44  »  *  »  *  J  ♦♦■»»  -f  *  4  4  *  4  4  »*•*  4  ■*  44  4  4  4  4-  4  4  -•  »*4M4  +  ¥<4*-»4-4444+V**¥**»-4****< 

Ic (EU.37.1) GO  TO  50 
GO  TO  500 

229  °RINT*,"M3  ELLIPTICAL  INTERCEPT  OR3ITS  POSSIBLE  FOR  GI7EN  DATA’* 

ST  0° 

50  I57  (KK.E0.1  )  GO  TO  111 

77  T  =  T  + . 1 

IF  (T.GT.20. ?G0  TO  229 
GO  TO  73 
5  C l  CONTINUE 


044»4*4"«4-4>44**-*:4*4-«4>*4 


»44*»4*4»  4  4  4  44*» +  ‘ 


4-  «•**■»  **** 


0 E 7 E L 0 3  the  STAtE  TRANSITION  MATRIX 


»  44  4  44»  4  4  4  4  »*  >44 


JY  =  1 
LL  =0 
I  - 1 

T  W=  ,  0 0 3  1 
XW=X+TW 
YW=Y+TN 
7W=7+TW 


7iw=vn*TW 

V2W=VI2+TW 

V3W=7I3-*-TW 

RI 2=R2 

RI3=R3 

VIS1=VT1 

71 S2=VT  2 

7IS3=VT  3 

XC  =  XW 

YC=Y 

ZC  =  7 

V1C=VI1 

V2C=VI2 

73C=VI3 

R0s$9Rt (XW«*2+Y*  *2  *-7**  2) 
VC=SORT (Vri*‘2+Vl2**2*VI3*'2> 
DOT  =  XW-  VI14-Y*7I  2  *2 4  VI  3 
A=l./ <70A*  2-2. /PO) 

A=A=»sm 

CALL  FINOYf  XX,  OPT, P.2,  VO, 7  O*  }> 
CALL  LOCATE (XX ,rO»VO»FtG»3»C2j 
Rl=c»XC+G' 7iC 
R2=F* YC+G1  V 20 
R?  =  F*7C*C’ 7  30 

R=SORT(RV  *2.«-c’?**2.*R3**2.) 
OF=XX* (7>» S-l . >/<R»RQ> 
D&=i.-XX«* 2,‘CC/R 

71  =  OFj*KC  +  DG*71C 

7 2  =  3F*  YC  +  C'G,F'/2C 
V3snp47r fnG*V3C 

7=  CORK  7t*  *2.+V2<‘*2.  +  V3*"-2.> 
CHFCK  (1,I)  =  (Rl-FID/TW 
CHECK  (2,1)  =  (22-t  I2)/TW 
CMECK<3,I)= ( ?3-rI3)/TW 
CHECKS ,I>=  (71-VIS1)  /TW 


CHECK  (5,1)  =  (V2-VIS2)  /TW 
CHECK  (S  ,  I)  =  (V  3- VI  S3)  /TW 
Ic  (LI  .£0.1i GO  Tr  96 
.IF  (I.  EQ.l)  CO  TO  7  2 
I c  ( 7  ,  £0 . 2)  GO  TO  73 
IF(I.E9.3)GO  TO  7k 
IF(.T,c9.4)  GO  TO  75 
IF (I . EO .5) SO  TO  76 
IF  (1, EO. 6) GO  TO  9b 
1=2 

R9=S0RT  (X*1*  2  +  YW"  +  2 +7  +  * 2) 

DOT=X« VI1+YW*VI2+7*VI3 

xc=x 

YC.  =  YW 

GO  TO  19 

1=3 

RO=SQRT  (7W*  *•  7+Yr  *  2+X»*  2) 
09T=X* VIi+Y»VI2+7W*VI3 
YC=Y 
7C  =  ?W 
GO  TO  99 
1  =  4 

RO  =  SORT  (X*»24-Y**2  +  7**2) 

VO = SORT  (VlWr‘2*VI2**  2+VH"’) 

D0T  =  X*71W  +  Y*-/I2*?*VI3 

ZC  =  7 

VlC=VlW 

Gu  TO  99 

1  =  5 

VC=SOFT  (VIi',*2  +  V2W*»2*VI3»»2) 
90T=X',VIH-Y*V2W  +  T-»V13 
V1C=VI1 
V2C  =  V2'( 

GO  TO  99 
1  =  6 

VC  =  SORT (VI 1**2  +  VI Z** 2+V3W*  *  2) 


AIR  FORCE  INST  OF  TECH  WRIGHT-PATTERSON  AFB  OH  SCHOO— ETC  F/6  16/2 
ERROR  PROPAGATION  FOR  6R0UND-T0-SATELLITE  ELLIPTICAL  INTERCEPT  — ETCCU) 
DEC  78  M  J  VON  PLINSKY 

AFIT/GA/AA/78D-8  NL 


OOU  TtWOGO 


00T=X^VI1+Y*/I2  +  '?»V?W 
V2C=VI2 
V  3C=V3R 

LL  =  1  i 

GO  TO  99 
9C  CONTINUE 

444  4.-44  4  4444  44  4  444444*444 

C  COMPUTE  UP3ATE0  COVARIANCE  MATRIX 

C*4*  *  44*4444444444*44*  44444444*44444444  4  4»  444V  *  4  4  4  4  4  4  4*4  444*  4  4  44  **4*4l 

CALL  MULT(CHE3K,C0V,C,6,S,a> 

DO  7  1=1, f 
00  7  J=l, 5 

7  9(1,  J)=CHECK<  J,  T) 

CALL  MULT(C,T,PWI,6,6,6) 

4  444**4444444*4 4#  4  4  4*  4*444444* * *  4  4  444*  4  4  »  44  *  •  *4444444 44 4444444*44 44 44 

EXTRACT  POSITION  SU®-ELEMESTS  OF  30V ARIANCE  MATRIX 

4444^44*44  4444?  4  4444*  4444444*4  4  4  4  4  4  444  4  4  f  »  44  *♦  4444444444444  4  444*4444^1 

DO  2  1=1,3 
DO  1  J= 1, 3 
RR  (I » J) =°HI (I , J) 

CONTINUE 

44*44 4*44444 44 « 44*»* >4* 444 444*4444444444»444*4 44 44 444444 44 444444444 V4 

ALIGN  COVARIANCE  ELEMENTS  ALONG  SATELLITE  TRACK 

4**4 44 44*4*44*»4***4 *444*4 444444444 444, »»*4*  4*4*44444 44  4*444 44  4  44441 

VXR=VTl-V(« 

VYR=VT2-VE 

V7R=VT3-VE 

VK2FQPT <VXR* *2  +  VYR**2+V7R*»2) 

VXY=SORT(VXR** 2+VYR** 2) 

CPX=VXR/VXY  ; 

SPX= VYR/VXX 
P2=AT  AN (V7R/VXY) 

RCT(1,1)=C0S(P7)4CPX 

R0T(1,2)=C0S(?7)*SPX 

R0T<1»T)=-SIM(P‘,’>  j 

ROT (2,1) =-5PX  ! 


it 

J 


/4-1 


o  o 


R0T(2,2)=:?X 
ROT (2,3)*n. 

RUT <?»1)=$IN(PZ)4CPX 
R0T(3,2)s£IN(P7)»SPX 
ROT (3,3)=rOS(°2) 

CALL  MJLR { ROT  »  PR, 03, 3 , 3, 3) 

00  8  1*1,3 
00  8  J*1 »  3 

8  ££  (I,  J>  =R3T  (1,1) 

CALL  KULR<00  »  EE,££,3,3,3> 

00  11  T=l,2 
00  10  J«l,2 

10  CCVRU,  J)  =  E£(I  +  1,  UH> 

11  CONTINUE 

COMPUTE  ERROR  ELLIPSOID  RESJuTS 


*  *  J  *  +  *.«..*******>!  »  ***** 


17=2 

CALL  El  GRP  (COUP.  ,2»2,2»W,‘t,IZ,WK*I£R> 

OCOsVA’VTl/  ('/*  VVS)  *V54VT  ( ’/*  V  VS)  M/£*  V TVIV+WS) 
MANr>*ACOS<300> 

HANG= HANG* 18 3 . / 3. 141592GE353 

L1  =  V4<1) 

L  2=P  ( 2) 

IMLi*L2..r  ,1.)  GO  T3  111 
LI “6378 1-Vi .  “SORT  (LI) 

L2=63731!»:t*>ORT(L2) 

TM=T' 13. '4^6334?? 
b6  CONTINUE 

CEP=S0RT(L1,'L2> 

OVI  =  Vl*7.9053f-828 
LK=LKH 

PR  I  JT4,  ” - **" 

?«rn*  #**TlMt  TO  INTERCEPT  =  ",M,*,MH>  *‘»T »”(T0) 

PRINT4  ,  "EL  LI3  5  01 D  AXES  LE'IIHS*  “  , -1  “ 

PRINT4, "  •*,L2,*'(M)  M 


PRINT*  ."INTERCEPT  ANGLE=",  H  AN6  ,  "  OE5 1  ” 

PRINT',  •'INTERCEPTOR  LAUNCH  'll  w  OCI  F  t- OVI  ,**( KM/SEC)  ** 
PRINT*, ***'  ‘C£P  =  ",CEP,“(M>" 

TF(L<.ST.r0)"O  to  27 
CC=CD*1 
I=2*CQ-1 
j=2*cn 
YA  (CD)  =  L1 
ZA (GO) =  L2 
TT(CO)=T 
ARY (I)=L1 
ARY(J)=L2 
OEP(CD>=CEP 
ICO  T=T ♦IF* .325 
KK*1 

JJJ=JJJ*1 
GO  TO  76 

27  CALL  PLOTR<YA,ZA,TT,ARY,C3> 

IF  (III, EO. 1)50  TO  66 
GO  TO  111 

66  CALL  ERROEP,TT,CO> 

GO  TO  111 

222  IF  (<LM.E0.12)G0  TO  223 
111=0 

PRINT* ,  **  " 

PRINT*,"  " 

PRINT*  ,  **  " 

GO  TO  111 

223  CALL  PLOTE 
STOP 

END 


SUPROUTTJE  L3TIME 
COMMON  THETA 

COMHON/LSTIHE/LONGI,DAYf  THE,  THETA GD 
REAL  LCNGI , LONG 
°I=3. 141592654 
LONG=LON5I*Pt/ieO. 

TIME=TIM£/2M0. 

DAY=DAY*TIME 

THETAsTHETA  GO  41. 0  0  27  375(P3»2.  *  PI'DAY+LONG 

RETURN 

END 


subroutine  srrEfx,Y,z,vx,i/y,vz) 

COMMON  THETA 
COMMON/SIT £E/LA*rI,  HI 
REAL  LA T *  L  A  TI 
AE=i.C 3 3005156 
HE=.9956"10553 
£=.39131 
ei=l. 1L 1592674 
LAT=LATI*PI/180. 

H=HI/2.0923572?TE7 

XI  =  (AE/ <1.-E**2.»SIN(LAT) *  *2.  I  *♦.  SH>*COS  CLAT) 

?=(AE»(l.-S**2.>/U.-E**2.*SIN(Ur>**2.>**.5+H>*SIN<LAT> 

X=XI*C3S(TH£TA) 

Y=XT*SIN (THET  A ) 

'/EARTH=.JfS833fC01 

VX=-VEARTH*Y 

VY=V£ARTH*X 

V7=0. 

RETURN 

END 
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SUBROUTINE  M  >LT  (0»B>C>L»M|'O 

DIMENSION  0(5, f.)  , 9 <6 ,6) ,C ( S, > )  , AO ( a, 6) , 90 <6 ,6) ,C0 (6 ,6) 
ICO  DC  10S  J=i, H 
00  104  1=1,  L 
1C4  AO(i, J) =0(1, J> 

DO  10£  K=1,N 
108  30  ( J, K)  =*«(J,<) 

00  112  1=1, L 
00  112  J*1,N 
CO (I  *  J) =9 • 

00  112  K=l,4 

112  CO(I, J)=CO(I, J) +AO(I, K>*0){<,  J) 

00  116  1=1,1 
00  116  J=l, N 
116  C(I*J>»C9(Zti> 

RETURN 

END 


I 


i 


r 

t* 


t 
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I  SUBROUTINE  H  it*.  (A,  P,0  ,L,M,  Ml 

DIMENSION  M1,3),o(3,3),0(3,3),A*(3,3),*C<3f3),C9<3,3) 
{  2UC  00  208  J*1,M 

,  00  204  '1*1,  L 

!  204  Ak(I,J>=UI,J> 

no  2 08  K*i,N 

208  3C  (  1,  K)  =P(J, <) 

00  212  1*1, L 

DO  212  J=1,N 

CB(l,J)=Ot 
00  212  K*1,N 

1  212  C"(I,  J)  *Or(I,  J)  ♦AP«,K)*8;<K,  J) 

'  <  00  216  1*1, L 

1  90  216  J=1 , N 

216  0(1, J)=CO(I,J) 

i  RETURN 

END 


i 

i 


1 

! 


i 
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SUBROUTINE  OR3IT(XI,YI,ZT,  ML,  VYI  ,\/ZI,RO,  VO, A) 

DIMENSION  0(5,6) ,B(6,b) ,0(5,5) ,0(3,1) ,SS(3,1> 

REAL  I 

RE  AO*,  A,  E,  I  ,OMEGA, WMEGA,AN3MY 
IF (EOF { yLI  NPJT ) .NE.O) STO" 

PRINT* ,  "SAT  £LL  I TE  *  S  ORBIT*.  ELEMENTS  ( A , E , I , OMEGA , WMEGA , ANOMALY > I 
PRINT*,  ••***«****•#  “»A,C, t ,  OMEGA,  rfl EGA ,  A NOMY , ** 

P=A*(1.-E‘  *2) 

Z=3.1G15>92653?R/180. 

I  =  I*Z 

OMEGA=OMEGA*Z 
WMEGA=WMESA*Z 
ANOMY=ANOMY*Z 
R=PM1.+E*  COS  ( APOMY)  ) 

B(l,l)=R*rDS(ANOMY) 

3(2,1)  =  R*FIN( A  NOMY) 

3 ( 3, 1)  =  0, 

D(l,l)=-(l./P**.5*SIN(AN0Mf)  ) 

D(2,1)=1./3**.S' (E*COS (AN31Y)  ) 

0(3,1)*0. 

3(1,1)=C0S(0MEGA)*C0S (WME35 ) -S IN ( DME GA ) *  SIN ( WMEGA) *COS ( I) 
a(i,2)=-C03(0MEF-A)*?IN(WME3A)  -  SIN  (  D^EGA )  *  COS  ( WMEGA )  *COS  (I) 
Q(1,3)=SIN'(0MEGA)*SIN(I) 

0  (  2, 1)  =  SIN(  OMEGA )  *COS  (WM-  31 )  OOS  (OMEGA)  -SIN  (WMEGA)  *COS(  I) 

0(2,2)  =  -SIN  OMEGA  )*SIN(  WMEGA)  ♦  COS  (  DM  EGA)  *  COS  ( WMEGA)  *COS(I) 
O(2,3)=-C0S()MEGA)*SIN(I> 

0(3, i)=SI P (WMEGA )* SI N(I> 

0(3, 2)= CO F( WMEGA)* SI N(I) 

0(3, 3)  =  COS  ( I) 

CALL  MJLT (0,B,C,3, 3,1) 

DO  11  <=1,3 
SS(«,1)=C(<,1) 

3(K, !)=')(*’,  1) 

CALL  M'JLT  (0,B,C,3,3,1) 

XT=SS(1,1) 

Y1=SS (2, 1)  . 


\ 


A-i+ 


?i=ss  n,t> 
vxi=c<i,i> 

VYI=C(2,1) 

V7I=C (3  *  1) 

RG=(XI**2.*YI**2.  +  71**2.)  +  ».5 
V0=(VXI*+2.  +YYIJ  2.+VZI*'2.) ►* .5 
A=i./(VO** 2-2./F0) 

A=A<3S<A) 

RETURN 

END 
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SUBROUTINE  FIMOX(XX»DQT»RD»VD  ,  IND»  A) 
COMHON/TIME/T 
X  =  T/A 
1  =  0 

1  Z=X>+2./A 

C=tl.-30S(S0*f  C7)>)./Z 

S= (SOFT  CM -SIN (SORT(Z)) >/S3RT (2**3.) 

TT=QOT*  X*  •  2  .*C  +  <l.-R0/A)~<»  +  3.  +  S+*:>*-X 

0T=X*-*2*C+30T*X*  ( l.-Z*S)  R D ♦  <  1 

3ELT=T-TT 

1  =  1+1 

I F  < I • El • ?  0 )  G3  TO  10 
IF(ABS(OE-T)  ,wT.l.E-6>G0  T3  20 
X=X+DELT/3T 
GO  TO  1 

10  PRINT*  » "CONVERGENCE  NOT  <-H  3E--X=" ,  X 

20  XX=X 
INO  =  I 
RETURN 
ENO 


A-t<* 


f 


SUBROUTINE  LDCATE(X<.RO»'/3,L',GfSt- 

COM-.ON/TIME/r 

7=  XX**2  . /£■ 

S=(SO^T  (T)-SrN(SQP.T(Z)))/53?r  12**  3 

C*(1.-30S(S-T»T<7)  )  )/7 

F=l.-XX**2*C'RO 

G=T-XX**3.*S 

CC=C 

ZZ*7 

RETURN 

END 


4 

4 


:,Z7,A> 

.  > 
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SUBROUTINE:  GAUSS (NU, MAGI,  -US?,  F,G,03,LM) 

CUMNON'/TI*F/T 
REAL  N'J ,  M  f  3 1 ,  M  A  0  2 

AA=SGRT  (Mt5 1*  HAG2)  *SI  N  (N'J)  /  30  R  T  Cl . -0  OS  (N*U>  ) 

LM  =  0 
Z  =  G. 

1  =  0 

1  I«r(7.GT..25)G0  TO  55 

IF (T.lT .-.25) GO  TO  57 
GO  TO  59 

5f  C=(i.iCOS(SORT  {"')))  77 

S= (SOFT (7) -SIN (FORT ( Z ) ) ) /SORT (7**3) 
DC=Ci.-Z*S-2.*C)/(2.*7> 

OS=  (C-3  .*S)  / (2 .*  Z) 

GO  TO  50 
57  C= (l.-OOSM(-T) )/Z 

S= (SINH (SORT ( -7) ) -SORT  <-T))/SQRT(<-7)**7) 

0C=l./24.+2  .*'T/~2G.-3.*Z*,2/4C  3  20  .  f  4 .  *7’-  *  3/36288  0  0  . 

□  S  =  l./120.+2.*7/5  04Q-3.*7*'  2/3  62630. +‘  .  2* *  3/399168  00 . 
GO  TO  50 

59  C=l./2.-7/24.*7'*2/?2  0.-,'»»3/4  0  320. 

S= 1  ./3. -7/1 20 . 4-  *"**2/5  04  0.-  3/332930  . 

OC  =  l»/24.+2.'*7/~2G.-3.*Z't*2/40  325 .* 4 .• Z* *  3/362  8  80 0 . 

OS  =  l./l  23.  *2.*7/5Q40.-3.*?»*  2/3623  3  3  . +u .* Z** 3/39916800. 
GO  TO  50 

5  c  i=in 

t,  Y  =  MAG1*MAG2-AA*  (l.-?*S)/S0RT(0) 

IF(Y.LT.O)GO  TO  58 
X=SORT(Y/C) 

TT=X*M*34AA*50rT  (Y) 

.  DELT=T-TT 
OELS=AOS( PELT) 

IF (0EL3.LT . i.E-C )GO  TO  10 

OT  =  X**3*C?S-7.*S*DC/<  2.*C>  )+AA/8.M3.*S'  SORT  (Y)  /C*AA/X) 

Z=Z+OELT/DT 

IF  (I.GT .500)30  TO  80 


.  ^-/8 
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GO  TO  1 
F=l.-Y/MA3i 
IFC7.LE.3.) GO  TO  88 
G=AA«SORT(Y) 

QG=1  •  -Y /MAG  2* 

LM=2 

GO  TO  33 
6  C  LM  =  0 

GO  TO  33 
86  LH=1 

83  RETURN 

ENO 


/*-/? 


SUBROUTINE  PLO^FfYAjZA^Tj^^CD) 

INTEGER  CD,C° 

01  PENSION  YA(40D)  ,26  (400)  »TT  (’4  0  0)  , A?Y(81Q) 

CDi=C04i 
CO  2*C0*2 
CP=2*C3 

CALL  PL0T(B  .,-2. 5, -3) 

call  plot  ( o  • ,  i  .  zz  £>,  -3 ) 

CALL  PLOT(0  .,3  .Ti>,-2) 

CALL  PL0T(6.2€,C.,-2) 

CALL  PLOT ( O . , -8  .75,-2) 

CALL  PL OT  ( - 6  •  ?1> ,  0  •  , - 2 ) 

CALL  PLOT (.75, •  r5  * -3) 

CALL  SCALF(TT,5.,CD,i) 

CALL  SCALE( ARY ,7 .  25»CP»1) 

YA (C041)=ARY(CP*1) 

YA (CD+2)=6RY(CP+2) 

ZA (CP41)=£RY(CP+1) 

7A  (Cr>42)=ARY(C°42) 

CALL  AXIS(0 0 . ,19HTIME-0'-RLlGHT  (TU) , -19 ,5 . , 0 . ,TT (COi ) , TT (C02) > 
CALL  AKIS(0.,Q., 21HCROSS  RACK  ERROR  CM ) , 21 ,7 . 2b , 90 . , ARY (CP *1 ) , ARY 
MCP*2)) 

CALL  LINE(TT,YA,CO,l, 0,0) 

CALL  LINE(TT,7A,CD,1, 0,0) 

RETURN  ,  ' 

ENO 


v 


i 
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SUBROUTINE  ERR(DE°,TT,CO) 

DIMENSION  DE°  ( A  PCI )  ,TT  (409) 

INTEGER  CD 

CALL  PLOTO.,-2.5»“3) 

CALL  PLOT(0.,  1.37b, -3) 

CALL  FLOT<0.,3.75,-2> 

CALL  PLOT(6.?b,C.,-2) 

CALL  PLOT(0 -8  .75,-2) 

CALL  Fi.OT(-5.25,0.,-2) 

CALL  PLOT  <  .  73,  .75  , -"*)  S 

CALL  $CALF(TT,P.,CD*1)  ] 

CALL  SCALE(OE»,?.25,CC,l) 

CALL  AX  IS  (  0  • ,  0.,  19HT  I ME-OF- FLIGHT  (  T  U)  , -19 ,5  .  ,  0  . ,  TT  (CD  +  l)  ,TT(C0*2| 
*> 

CALL  AX  ISO. ,0.,  34HCIRCUH  RUED  CROSS-TRACK  ERROR  ( M)  ,  34, 7 . 25, 3  0. 
J  DED  (CD«-1)  ,0Eo(Cr  +  2)) 

CALL  LINE<TT,OE°,CO,1,0,0» 

RETURN 

EHC 
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